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Abstract. We study experimentally phase separation and transient gelation in a 
mixture consisting of polydisperse colloids (polydispersity 6%) and non-adsorbing 
polymers, where the ratio of the average size of the polymer to that of the colloid 
is ~ 0.062. Unlike what has been reported previously for mixtures with somewhat 
lower colloid polydispersity (« 5%), the addition of polymers does not expand the 
fluid-solid coexistence region. Instead, we find a region of fluid-solid coexistence 
which has an approximately constant width but an unexpected re-entrant shape. We 
detect the presence of a metastable gas-liquid binodal, which gives rise to two-stepped 
crystallization kinetics that can be rationalized as the effect of fractionation. Finally, 
we find that the separation into multiple coexisting solid phases at high colloid volume 
fractions predicted by equilibrium statistical mechanics is kinetically suppressed before 
the system reaches dynamical arrest. 
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1. Introduction 

Well-characterized colloidal suspensions are often used as experimental 'models' of 
phenomena and processes occurring in condensed atomic or molecular systems such 
as crystallization, various liquid crystalline transitions and glassiness [HE]- While this 
'knowledge transfer' from colloids to the atomic or molecular world works well in many 
instances, there are also phenomena that are unique to colloids. In particular, real- 
life synthetic colloids are inevitably polydisperse: the suspended particles are never 
identical, but have a distribution in properties such as size, shape, charge, etc. Even if 
many properties are uniform, e.g. in a suspension of uncharged spheres made of a single 
material, there is an inevitable polydispersity in size. We define this polydispersity, a, as 
the standard deviation of the radius distribution divided by the mean. A bulk sample 
of polydisperse colloids is, for all practical purposes, an 'infinite component system'. 
In contrast, atomic or molecular materials almost always contain a finite number of 
components, and the species within each component are (apart from unimportant 
isotopic variations) identical. 

In colloidal experiments aimed at elucidating generic condensed matter phenomena, 
polydispersity is usually treated as a problem to be minimized. In such research, the 
aim is to use particles with as low polydispersity as possible. For example, micron-sized 
fluorescent silica particles for confocal microscopy can be synthesized with a < 2% [3]. 
A low degree of polydispersity does not seem to introduce qualitatively new phenomena 
in many instances. Thus, experimentally [4], slightly polydisperse hard-sphere colloids 
(er < 5%) show the equilibrium phase behaviour predicted by computer simulations for 
monodisperse hard spheres [5], viz., crystallization at volume fractions (ft ^ 0.5. Even 
in this simple case, however, polydispersity can have significant quantitative effects. 
For example, theory predicts [6l [7] that even at a as low as 3 or 4%, the melting and 
freezing volume fractions (4>m and (ftp) may be shifted to higher values compared to their 
monodisperse values of 0.545 and 0.494 respectively. While such predictions are non- 
trivial to verify because of the difficulty in measuring volume fractions in polydisperse 
samples, experiments unambiguously demonstrate [8] that polydispersity significantly 
affects kinetics: nucleation was delayed dramatically on increasing a from 4% to 6.8%. 

Predicting the phase behaviour of polydisperse systems in general, and of 
polydisperse colloids in particular, poses a significant theoretical challenge. In few- 
component molecular systems, equating the chemical potentials and pressures of 
coexisting phases is a practical method for calculations. But such a 'brute force' 
approach becomes increasingly cumbersome as the number of components grows, and 
entirely new approaches are needed to treat properly the polydisperse, or infinite- 
component, limit. Significant advances have been made in recent years (reviewed in 
P,[in]), especially on the 'reference system' of polydisperse hard spheres. 

In monodisperse hard spheres, increasing (ft of a fluid leads first to fluid-solid (where 
'solid' means crystal) (FS) coexistence in the interval 0.494 < (ft < 0.545, and then the 
equilibrium state is single-phase solid from (ft = 0.545 to close packing ((ft m 0.74) j5]. We 
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denote this sequence in the self-explanatory notation as 'FS — > S'. Theory predicts [7] 
that a small degree of polydispersity changes the sequence to FS —> S — >■ SS, i.e. at high 
enough 0, we should find fractionation into coexisting solid phases with different lattice 
parameters. Above cr ~ 6.5%, the single phase solid (S) disappears, and increasing 
beyond the fluid-solid coexistence region brings the sequence 'FS — > FSS — > SS —> SSS'. 

Theoretical predictions have also been made for the effect of polydispersity on the 
equilibrium phase behaviour of mixtures of hard spheres and non-adsorbing polymers. 
In such a colloid-polymer mixture, the exclusion of polymer between the surfaces of two 
nearby particles gives rise to an effective 'depletion' attraction whose range and depth 
are proportional to the polymer size and activity (or, equivalently, chemical potential) 
respectively [llj . Since either the polymer or the colloid or both may be polydisperse, 
a rich variety of effects may be possible in principle. 

The equilibrium phase behaviour of mixtures of monodisperse spherical colloids and 
polymers is well known (reviewed in [12j[T3]) from both theory [T4"l[T5] and experiments 
[El [T7j, and depends on the ratio of the colloid radius (a) to the size of the polymer as 
measured, e.g., by its radius of gyration (r g ), £ = r g /a. At low enough £, the region of 
fluid-solid coexistence found at (f>p = 0.494 < < (pM = 0.545 for monodisperse hard 
spheres broadens out. At high enough £, however, the phase diagram additionally shows 
regions of gas-liquid and gas-liquid-solid coexistence. The cross over occurs at £ = £ c , 
where theory predicts £ c ~ 0.3 [HI [13] and experiments find £ c ~ 0.25 [IB] . 

The effect of polymer polydispersity in such mixtures has been studied theoretically 
[18] . Calculations uncovered no new qualitative effects. In the rest of this paper, we 
therefore do not consider further the effect of polymer polydispersity. 

Theory suggests that the effect of colloid polydispersity in colloid-polymer mixtures 
is richer. Fasolo and Sollich [19] have published schematic representations of the effect 
of increasing colloid polydispersity, a, on colloid-polymer mixture phase diagrams at a 
range of £, but have made detailed calculations only for £ > £ c . Experimentally, limited 
data are available for mixtures with (£, cr)= (0.4, 0.18) [20] and (1, 0.54) [21] . In both 
cases, significant colloid fractionation was observed in coexisting colloidal gas and liquid 
phases, with the degree of fractionation obeying a universal law [20J . 

While (of necessity) all previous studies of colloid-polymer mixtures in the £ < £ c 
regime have used polydisperse particles, we are aware of no attempt to date specifically 
to study the effect of colloid polydispersity on the phase diagram at these smaller size 
ratios. In this paper, we report the first results of such a study using a hard-sphere- 
like suspension with o > 6%. We find agreement with certain aspects of the phase 
diagram topology predicted by Fasolo and Sollich [19] , but that separation into multiple 
coexisting solid phases is suppressed. These results lead to a slight revision in the 
interpretation of previously published data on multiple glassy states in a colloid-polymer 
mixture with a somewhat higher a [221 123] . We also detect the presence of a metastable 
gas-liquid binodal, which gives rise to multi-stepped phase transition kinetics. 
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2. Experimental materials and methods 

We used particles of polymethylmethacrylate (PMMA) sterically stabilised by 
chemically-grafted poly-12-hydroxystearic acid (PHSA) |24J synthesized in house and 
suspended in a mixture of cis- and trans-decalin (density 0.88123 g/cm 3 ). We 
characterized these particles using small-angle X-ray scattering (SAXS) (reviewed in 
[25J). SAXS measurements were carried out at the ID02 beamline at the European 
Synchrotron Radiation Facility in Grenoble (France). Two sample-to-detector distances 
of 10 m and 1 m were used to cover a scattering vector (q) range of 0.008 nm -1 < 
q < 3 mil" 1 . Colloidal suspensions were contained in a flow-through capillary cell of 
diameter 1.8 mm. Measured two-dimensional scattering patterns were normalized to 
an absolute intensity scale and then azimuthally averaged to obtain the intensity as a 
function of q, I(q) [26J. Fig. [TJ shows the background-subtracted SAXS intensity from 
a dilute suspension. The data can be adequately described by a polydisperse core-shell 
model with a Schulz size distribution [26]. The best fit parameter values are: mean core 
radius and polydispersity 152 ± 2 nm and 6.1 ± 2% respectively, and volume fraction 
0.006. 

The fitted mean thickness of the PHSA hair was 6 ± 1 nm, agreeing with the value 
obtained by Cebula et al. from core-shell fitting of small-angle neutron scattering data 
|27J and viscosity measurements (28]. Our core-shell fitting (and that in reference |27j ) 
assumed a homogeneous shell, while viscosity measurements return a value relevant for 
hydrodynamic properties. Cebula et al. also measured structure factors of concentrated 
suspensions, and from the fitted effective hard-sphere diameter deduced a shell thickness 
of 10 nm. They attributed the difference with core-shell fitting the solvation of the 
shell and the polydispersity of the PHSA chains. Indeed, chains as long as 15-20 nm 
may be present [28]. The longer chains will have a significant effect in determining an 
effective hard-sphere diameter, which presumably is the relevant characterization for 
our purposes. We therefore take the shell thickness to be 10, with a polydispersity of 
±4 nm. This returns a final estimated particle size of 162 nm and a polydispersity of 
6.2% (the latter figure having been obtained by adding the variances of the core and 
shell). 

The polymer was linear polystyrene (PS) with mass-averaged molecular weight 
M w = 114,200 daltons (Aldrich GPC Standard, unknown polydispersity), with an 
estimated radius of gyration of r g 10 nm at room temperature, where we expect 
PS to behave as nearly-ideal coils in our solvent [29J. Thus the (average) size ratio of 
our system is £ ~ 10/158 = 0.062. We estimate that the overlap concentration, c*, is 
reached at a polymer concentration of c p ~ 45 mg/cm 3 . The highest concentration we 
used was < 0.4c*. 

We calibrated the colloid volume fraction using the theory of Fasolo and Sollich for 
the phase behaviour of polydisperse hard spheres [7]. The latter predicts that in the 
vicinity of a ~ 6%, the solid branch of the FS coexistence region displays a maximum 
at (j) = 0.58. We assigned (f> = 0.58 to a pure colloid sample with close to 100% crystals 
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Figure 1. Form factor of PMMA colloids (points) from SAXS fitted with a 
polydisperse core-shell model convoluted with an instrument resolution function of 
width 0.0026 nm" 1 (continuous line, red online). 



(■ on the c p = axis, Fig. [2]). This procedure introduces systematic uncertainties into 
our values in the region of 1% or more, although the statistical uncertainties in the 
quoted values are considerably smaller. 

Mixtures were prepared in 1.8 ml cuvettes by adding solid polymer to colloidal 
suspensions with known and tumbling overnight. Phase behaviour was studied by 
direct visual observation over 3 to 14 days at 25°C. Time-lapsed images were also taken 
in a number of cases. After observation, samples were homogenized again by tumbling, 
and then diluted to yield the next generation of samples for observation. All our data 
points therefore lie along 'dilution lines' in the (0, c p ) plane. For convenience, however, 
we describe our observations below in terms of the progressive addition of polymer to 
samples. 

3. Results 

The different kinds of behaviour we observed are summarised in a phase diagram in 
Fig. [2l At < 0.55 and low enough polymer concentration, c p , samples remained 
single-phase fluids (♦ in Fig. [2]). Adding sufficient polymer brought about fluid-solid 
(FS) coexistence (□ in Fig. [2]), with iridescent crystallites nucleating throughout a 
sample and then falling to the bottom of the cuvette. Further addition of polymer 
gave rise to samples that separated into two coexisting, non-crystalline phases with 
different turbidities (and therefore volume fractions) in which the interface between 
the phases was sharp and remained horizontal when a sample was tilted (• and * in 
Fig. EJ). This is the behaviour expected of samples undergoing gas-liquid (GL) phase 
separation. A number of samples initially showing GL phase separation subsequently 
nucleated crystals from the denser, liquid phase, giving rise to final states displaying 
gas, liquid and crystal (GLS) phases (* in Fig. [2]). 
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Figure 2. Observed phase diagram, with legend given in the inset, where G, L and 
S stand for gas, liquid and solid (= crystal) phases. Three kinds of FS coexistence 
samples are distinguished, those with < 20%, > 80% or between 20-80% of crystals: 
S < 20%, S > 80% and 20% < S < 80%. Lines are guides to the eye of where various 
phase boundaries lie as deduced from the data points; the exception is the dotted line, 
which is a schematic of the attractive glass line. (Colour online.) 

At even higher polymer concentrations in this region, the behaviour was observed 
to change again (A in Fig. [2]). The final state of these samples looked at first sight very 
similar to those that had undergone GL phase separation. But the interface between 
upper and lower phases was often visibly rough, and did not tilt when a sample was 
tilted. In some cases, we obtained sedimentation profile from time-lapse images; these 
showed a delay time before rapid sedimentation began. Such 'delayed sedimentation', 
giving rise to a solid-like sediment, is characteristic of transient gelation [30] . 

In the region 0.55 < < 0.58 and c p = 0, we again found FS coexistence. 
The amount of crystals grew as <j) increased (EB), until we reached a sample that was 
essentially fully crystalline (■ in Fig. EJ), which delimits the high-concentration side of 
the FS coexistence region. Upon addition of polymer, the width of the FS coexistence 
region stayed roughly constant, but displayed a re-entrant shape. Interestingly, all 
samples to the immediate right hand side of this region failed to crystallise (0 and V in 
Fig. [2]). This is in striking contrast to what was observed in a colloid-polymer mixture 
containing colloids with a somewhat lower a of 5% [31] , where immediately to the right 
of the FS coexistence region, a region of full crystallisation was found. On the other 
hand, our observations are consistent with previous work on another mixture with a 
£ similar to the present one and a ~ 7% [22| 123] . In that work, it was assumed that 
all non-crystallizing samples to the right of the FS coexistence region were glassy [22] , 
although it was noted that this assumption might not be correct [23] . An unambiguous 
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Figure 3. Topology of the phase diagram of a colloid-polymer mixture with £ 
appropriate to our experimental system according to [19]. (a) a < 0.6. The bold 
parts show what is observed experimentally. Inset: a schematic representation of the 
phase diagram for a mixture of non-adsorbing polymers and monodisperse colloids, 
(b) 0.6 < a < 0.75. (c) a > 0.75. In all cases, we have added a metastable gas-liquid 
binodal (dashed curve). 

method for deciding whether a non- crystallizing sample is a fluid or a glass is dynamic 
light scattering, which measures the intermediate scattering function (ISF). The ISF 
decays completely in fluids, but has a finite value in the long-time limit in glasses. 
But our samples are turbid, and would require the use of specialised light scattering 
equipment to probe their ISF |32j. We therefore chose a less precise, but much easier 
to implement, visual method. Samples whose menisci stayed horizontal when vials were 
tilted were classified as fluids, while samples whose menisci tilted immediately with the 
vials and stayed tilted for some time were classified as glasses. 

Since we classified samples by visual observations alone, the boundaries between 
the different types of behaviour described above are not entirely sharp. Nevertheless, 
taking into account uncertainties in observation and in sample concentration, we can 
locate these boundaries to within about one data point in a typical dilution sequence. 
The approximate position of the boundaries demarcating the different types of behaviour 
are sketched in Fig. [2j 

4. Discussion 

We interpret our data within the framework of the possible phase diagram topologies 
in polydisperse colloid + polymer mixtures given by Fasolo and Sollich [19]. For low 
enough £, and certainly for £ < 0.2 (compare Figs. 2 and 5 in [19]) the topologies for 
low, medium and high polydispersities are shown in Fig. [3J (redrawn from Fig. 5 in 
[19]). From the sequence of phases expected for pure hard sphere (i.e. c p = 0) as a 
increases (see in particular Fig. 3 in [IE]), we may deduce that Fig. [3ta) should apply 
for < a < a cl « 6 + %, Fig. 0(b) for a cl < a < a c2 « 8~%, and Fig. [3(c) for a > a c2 . 
In other words, Fasolo and Sollich's theory predicts sharp transitions between different 
phase diagram topologies at well defined values of the polydispersity, a c i, a C 2, etc. In the 
vicinity of these 'critical polydispersities', a small change in a can give rise to qualitative 
changes in the observed phase diagram. 
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Note now from Fig. |^a) that when the polydispersity is low enough, o < er cl , the 
low c p portion of the phase diagram, picked out by bold lines in the figure, is topologically 
identical to that predicted for a mixture in which the colloids are monodisperse (inset, 
Fig. Efa)). Previous work on colloid-polymer mixtures at low £ and low a [T6J, I3T] 
confirmed this prediction. Addition of polymer was found to broaden the FS coexistence 
region. At higher polymer concentration, various non-equilibrium phenomena were 
encountered, which could be interpreted as due to the presence of a metastable gas- 
liquid (GL) binodal (dashed curve, Fig. |3]^a) and inset) 'buried' within the equilibrium 
FS coexistence region. 

No such broadening of the FS coexistence region is seen in our phase diagram, 
Fig. [2j Instead, the FS coexistence region remains more or less constant in width, 
although it has an interesting re-entrant shape. The behavior is consistent with the 
topology shown in either Fig. [3]^b) or (c). Since our measured polydispersity is just over 
6%, it is probable that Fig. [3]^b) is applicable. 

But whether Fig. E^b) or (c) applies in our case is in fact irrelevant for our purposes. 
In either case, the predicted equilibrium phase behaviour to the high-0 side of the FS 
coexistence region is not single-phase solid, but the coexistence of multiple solid phases 
with fluid. In recent simulations of polydisperse hard spheres [33], it was found that the 
locus of points beyond which crystallization failed to occur closely tracked Fasolo and 
Sollich's theoretical boundary [7\ separating the occurrence of a single solid (S) and two 
solid coexistence (SS) in the (0, a) plane. Apparently, fractionation into multiple solid 
phases is kinetically unattainable (cf. [8J). Our observations suggest that the same may 
be true in colloid-polymer mixtures: all samples immediately to the high-0 side of our 
FS coexistence region remained in the fluid state rather than show FSS coexistence. 

At even higher volume fractions, samples became glassy. The boundary separating 
fluid from glass appear to have a slight positive slope, suggesting re-entrance, i.e. the 
addition of polymer (= interparticle attraction) to the hard-sphere glass initially melts 
it into a metastable fluid, as found before J22J, 123] . 

At cf) < 0.5 we find a FS cloud curve (i.e. the low-0 limit of the FS coexistence 
region) with negative slope, consistent with any of the topologies in Fig. [3] and with 
previous experiments at similar |22j or lower [31J polydispersities. Inside the FS 
coexistence region, we find evidence for a metastable GL binodal, again consistent 
with previous work [161 122]. Such a binodal is known to be 'buried' within the FS 
coexistence region in the monodisperse system (inset, Fig. EJ^a)), becoming stable in 
its own right at £ > £ c [14] . A metastable GL binodal clearly continues to exist in 
polydisperse systems with low £. Inside such a metastable binodal (dotted lines in 
Fig. [3]), a homogeneous system may lower its free energy somewhat by separating into 
coexisting gas and liquid phases, although the absolute minimum in free energy (i.e. 
thermodynamic equilibrium) is only reached by separating into one or more of the solid 
phases indicated by the equilibrium phase boundaries (full lines in Fig. [3]). (For further 
discussion of 'metastable phase boundaries' of this kind, see [34J.) 

Compared to a system with lower a [31], the region of FS coexistence between the 
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FS cloud curve and the metastable GL binodal in our higher-cr system is very narrow 
(compare also [22] ). This is consistent with the finding that polydispersity destabilises 
the single-phase fluid relative to GL phase separation |35j , so that less polymer is needed 
to cause phase separation [19] . 

It is by now established that arrested phase separation within a GL binodal gives 
rise to gelation |36j. In colloid-polymer mixture with £ < £ c , such arrest is due to the 
intersection of the attractive glass line [22] with the GL binodal [36] , so that the gelation 
boundary is a GL tie line. In systems with low enough a, this intersection occurs very 
close to the GL critical point [36J, so that a region of complete GL phase separation is 
not observed [3U[36]. However, with the movement of the GL critical point to lower c p 
in a system with higher a, a region of complete (i.e. non-arrested) GL phase separation 
can be seen, as is the case in our system here. Significantly, this region is separated 
from the region of transient gelation (= arrested GL phase separation) by a straight 
line, Fig. |2) this is a tie line within the GL binodal. The high-0 end of this tie line 
indeed lies where we may expect the attractive glass line (shown schematically as the 
dotted line in Fig. [2]) to intersect the GL binodal . 

A number of samples just across the metastable GL binodal displayed interesting 
multi-stepped kinetics (* in Fig. [2]): they first separated into coexisting gas and 
liquid phases; subsequently, crystals nucleated out of the liquid phase, giving rise to 
samples displaying gas, liquid and solid phases. We interpret this as a non-equilibrium 
fractionation effect. The liquid phase is expected to have lower polydispersity than the 
parent metastable fluid [191 EQl [35] . This lower polydispersity is favourable for crystal 
nucleation. However, equilibrium FS coexistence cannot be attained, because a non- 
equilibrium particle size distribution was 'locked into' the liquid phase, and then into the 
crystals that nucleated from this phase. The (experimental) final state of these samples 
therefore contains out-of-equilibrium gas, liquid and crystal phases. These samples must 
be carefully distinguished from those showing equilibrium gas, liquid, solid coexistence 
in colloid-polymer mixtures with £ > £ c |16j . 

However, this two-stepped crystallization mechanism is only observed for samples 
close to the GL binodal. Here, GL phase separation is at its slowest, so that a higher 
level of fractionation should be kinetically attainable; this is favourable to crystallization 
from the liquid phase. Deeper into the GL binodal, GL phase separation occurs faster, 
so that we may expect a lower degree of fractionation on kinetic grounds. We suggest 
that this kinetic effect eventually renders the liquid phase not sufficiently fractionated 
to give rise to crystallization. Metastable GL coexistence is therefore the (experimental) 
final state of these samples. 

5. Summary and outlook 

We have studied experimentally the phase behaviour of a mixture of hard-sphere-like 
colloids and non-adsorbing polymers in which the colloids have a polydisersity of a > 6% 
and the polymer: colloid size ratio is £ ~ 0.062. The observed phase diagram displays 



Polydisperse Colloid-polymer mixtures 



10 




a number of features not found in similar mixtures with lower colloid polydispersityjjj 
Most notably, the FS coexistence region does not expand upon the addition of polymer, 
but rather takes a striking re-entrant shape of constant width. The physical origins of 
this re-entrant FS coexistence region is unclear at present. 

Inside a metastable GL binodal, we observed two-stepped crystallization kinetics: 
samples first separated into coexisting gas and liquid phases, and crystals subsequently 
nucleated from the liquid phase. We have interpreted this as a kinetic fractionation 
effect: a sufficiently fractionated liquid phase can nucleate crystals. Experiments in 
which the degree of fractionation in coexisting phases is measured, e.g. by fitting form 
factors from SAXS, can be used to confirm or refute this interpretation. Note that the 
multi-stepped kinetics we observed is consistent with a phenomenological model used 
previously to interpret observations from a colloid-polymer mixture with £ > £ c so that 
equilibrium GLS coexistence was possible [37J EE]- Moreover, if our interpretation of 
these observations in terms of fractionation is correct, then our results can be seen as a 
novel application of the Ostwald 'Rule of Stages' [39] in polydisperse systems. 




The constant-width FS coexistence region we observed is consistent with predictions 
of the phase diagram topology by Fasolo and Sollilch [19] for colloid-polymer mixtures 
with our £ and a. On the other hand, the regions of multiple solid coexistence predicted 
by this equilibrium theory are replaced in our experiments by metastable fluids or glasses 
- fractionated solid-solid coexistence is apparently kinetically inaccessible, either in 
mixtures of hard spheres and polymers (this work) or in pure hard spheres [33]. Such 
kinetic effects have not been predicted by any first-principle theory to date. Note that 
interestingly, coexisting hexagonal columnar crystal phases have recently been observed 
in mixtures of polydisperse colloidal platelets [40], although only after very long aging 
times (2 to 5 years). It is therefore possible that future experiments in hard spheres 
or hard spheres + polymer mixtures over similarly long time scales (years) will observe 
FSS or SS coexistence experimentally. 

Our intention here has not been to study the effects of polydispersity on the re- 
entrant glass transition known to occur in colloid-polymer mixtures with low £ and a 
[22], [23] . Indeed, such effects have not yet been probed in any detail, although one 
experimental report in a system with a ~ 23% indicates that such a high polydispersity 
may obliterate the re-entrance all together [H]. Our observations suggest that the re- 
entrance is still present in a system with a > 6%, although more data points and a 
detailed study using dynamic light scattering will be necessary to confirm this finding. 
It would also be interesting to probe the dynamics of the metastable fluids (0 in Fig. [2]) 
that seem to be prevented by kinetics from reaching equilibrium states in which multiple 
solid phases should coexist [i~9] . 
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% Note that the refractive index mismatch between particles and solvent here is closely similar to that 
in [16l|3T], so that differences are not due to differing residual van der Waals interactions. 
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